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Two Kerr-squeezed optical beams can be combined in a beamsplitter to produce non-Gaussian 
continuous-variable entangled states. We characterize the non-Gaussian nature of the output by 
calculating the third-order cumulant of quadrature variables, and predict the level of entanglement 
that could be generated by evaluating the Duan-Simon and Reid Einstein-Podolsky-Rosen criteria. 
These states have the advantage over Gaussian states and non-Gaussian measurement schemes in 
that the well known, efficient and proven technology of homodyne detection may be used for their 
characterisation. A physical demonstration maintaining the important features of the model could 
be realised using two optical fibres, beamsplitters, and homodyne detection. 
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INTRODUCTION 

Continuous variable (CV) systems provide flexible and 
powerful means for implementing quantum information 
schemes [JJ, in large part because there are mature and 
precise techniques for measuring the quadratures of light, 
most of which are familiar from classical communications 
technologies. Despite the need to deal with transmission 
losses, recent work has demonstrated useful distances 
comparable to those achieved with discrete-variable sys- 
tems [2]. Furthermore, research and development has 
progressed to the stage where CV quantum key distri- 
butions systems have advantages over discrete-variable 
methods [5H5]- 

The remaining stumbling block to the wider use of CV 
systems is that the most readily available CV systems 
and the the most developed detection techniques pro- 
duce only Gaussian statistics. This limitation rules out 
tasks such as entanglement distillation [BJ, quantum er- 
ror correction [7], and quantum computation. One way of 
introducing non-Gaussian statistics is through nonlinear 
measurements [8], but this approach negates one of the 
main advantages of CV systems, namely, the highly de- 
veloped technology that is available for performing Gaus- 
sian homodyne measurements. 

In this paper, we proceed along an alternative ap- 
proach, namely to use CV sources that produce non- 
Gaussian outputs. The importance of this area of re- 
search was shown recently by Ohliger et al. [5], who 
demonstrated there are serious limitations to the use 
of Gaussian states for quantum information tasks which 
may be avoided by developing useful and relatively sim- 
ple non-Gaussian sources. 

Non- Gaussian light can be produced by means of a 
%( 3 ) nonlinear medium, such as a single-mode optical fi- 
bre. Pairs of such Kerr squeezed beams can be combined 
on a beam splitter to produce CV entangled states, as has 
been experimentally produced using polarisation squeez- 
ing [10]. However, the non-Gaussian character of these 
states has not to our knowledge been explicitly demon- 



strated. 

In this paper, we use a single-mode anharmonic oscilla- 
tor |llj to determine the non-Gaussian entanglement that 
can in principle be achieved with Kerr-squeezed states. 
We characterise the non-Gaussian statistics through 
higher-order cumulants and gauge the level of entan- 
glement by calculating the Duan-Simon and Einstein- 
Podolsky-Rosen correlations. 

TESTING FOR NON-GAUSSIAN STATISTICS 

A Gaussian state can be most simply defined as a state 
with a Gaussian Wigner function, i.e. a state whose 
marginal distributions are Gaussian. For a CV state, 
the departures from non-Gaussian behaviour can thus be 
characterised by the skewness of the distributions of its 
quadrature moments, as revealed in nonzero higher-order 
cumulants [T2"] . 

We define the generalised quadrature X(9) at quadra- 
ture at angle 9 as 

X{9) = ae~ i0 + ah ie , (1) 

so that the canonical X quadrature is found at 6 = 0, 
with conjugate Y = X(^). 

For a Gaussian distribution, all cumulants higher than 
second order vanish, and therefore we can test for non- 
Gaussian statistics by a the presence of a nonzero third- 
order cumulant: 

k 3 (9) = (X"(9)) + 2{X{9)f - 3(X(9))(X 2 (9)). (2) 

While K3 7^ is a sufficient condition for non-Gaussian 
statistics, it is not a necessary one. In particular, k 3 
will vanish for a symmetric distribution in phase space. 
In the presence of such symmetry, the fourth-order mo- 
ment Ki provides the lowest-order test for non-Gaussian 
behaviour: 

k a {9) = (X 4 (9)) + 2(X(9)) A - 3(X 2 (9)) 2 - (X(9))k 3 (9). 

(3) 
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The fourth-order cumulant can be used to infer the 
negativity of the Wigner function [13) . which is con- 
sidered to be a direct measure of the nonclassicality of 
a state. It also allows comparison to the nonclassical 
states that have been experimentally demonstrated to 
be non-Gaussian, such as the number state [5] and the 
photon-subtracted squeezed vacuum [HI [15]. For both 
of these states, K4 scales quadratically with number. For 
the number-state for example, 



K4 = — 6n(n + 1). 



(4) 



In the analysis below, we will determine the regimes in 
which the Kerr-squeezed state is skewed to a similar level. 



NON-GAUSSIAN STATISTICS IN THE 
KERR-SQUEEZED STATE 

The Hamiltonian for the single-mode model, ignoring 
any effects due to loss and excess noise, is 

, 2 



H = h\ (a* a) 



(->) 



where \ represents the third-order nonlinearity of the 
medium and a is the bosonic annihilation operator for 
the electromagnetic field mode. 

For an input Glauber-Sudarshan coherent state, 
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where |n) represents a Fock state of fixed number, we may 
find analytical expressions for all the operator moments 
necessary to calculate the first four cumulants. 

The Heisenberg equation of motion for a can formally 
be solved to give 

-ixt(aa+a+i) a ^ j j 



a(t) = e 

whose expectation value in a coherent state is 
(o(i)) = 



ix£^Ja| (cos 2\t — i sin 2xt~ 1) 



(7) 



(8) 



Defining ag = ae l9 , we can write the mean of a 
quadrature moment as 

(X(6,t)) = (a e (t) + al(t)), (9) 

for which we already have a solution. The second moment 
is 



(X 2 (0,t)) = (1 + 20^+0^ + 0?) 



^ 2 



(10) 



where we have dropped the time argument on the RHS 
for simplicity. The third- and fourth-order moments are 

(X 3 (9)) = (4 3 + 3al 2 a e + 3ala 2 g +al + 3al + 3a g ) 7 

(x^e)) = (4 + 444 

+ Qal + 12a] 



6ai a e + 4ai a g 



6a 



•3). 
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To analytically calculate these moments, we use the fol- 
lowing expectation values and their complex conjugates: 



(4W) 


= a 2 e- 2U 


} e 


— 4:ixt A |c*| 2 (cos A\t—i sin 4x^—1) 
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= a 3 e~ 3U 


} e 


— 9ix*pl a | 2 ( cos i s i n 1) 


(44W) 
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-id 


e -3ix* e M 2 (cos 2 X t-i sin 2 X t-l) ^ 


(4(t)) 


= a 4 e- 4 *< 


>e 


— il§\t \ot\ 2 (cos 8xt—i sin 8x^ — 1) 
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i8x^g|a| 2 (cos 4xt— i sin 4%)t— 1) 
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These equations reveal several kinds of contributions 
to the dynamics, with different time scales. The sin and 
cos terms in the exponents can each be expanded, and 
for sufficiently small interaction time we can keep the 
first two terms in each, i.e. up to fourth order in time. 
We are left with a number of different contributions to 
the exponent. 

First, there is the nonlinear phase factor proportional 
to N\t, where N = \ct\ 2 . This mean-field frequency shift 
can be removed by a switch to a rotating frame, i.e. set- 
ting 6 = 6 + 2N X t. 

Second, the real exponent proportional to N)( 2 t 2 is 
responsible for squeezing, although in order to obtain 
quantum squeezing, i.e. below the coherent-state level, 
we also require the zero-point phase factors. 

Finally, there are the third and fourth order terms 
Nx 3 t 3 and Nx 4 t 4 , which for large N give the leading 
order contribution to the third and fourth order cumu- 
lants and hence are responsible for most of the skewness 
we see in the quadrature statistics. 

In typical Kerr squeezing experiments, the number of 
photons is large N ^> 1 in order to compensate for a 
weak nonlinearity x < 1 In this limit, we can derive a 
simple expression for the third-order cumulant of the Y 
quadrature rotating at the mean-field frequency, which is 
where skewness is most evident: 



K3 
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(13) 



The validity of this expression is demonstrated in 
Fig. [l] which plots the third-order cumulant for various 
photon numbers. The exact results for TV > 10 6 are in- 
distinguishable on this time scale from the simple cubic 



growth described by Eq. (13). Time is here scaled by 



N\ in order compare results that give the same Kerr ef- 
fect. On this scale, the third-order cumulant decreases in 
proportion to the square root of the number of photons. 
Note however, that the absolute size of the third-order 
cumulant increases with particle number, at a rate faster 
than (X) 3 : 



k 3 
(X) 3 



N. 



(14) 
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FIG. 1: (colour online) The third order cumulant, K3, of 
Y — X(ir/2) in a rotating frame as a function of time, for 
various photon numbers ranging from 100 to 10 6 , as labeled. 
In this and subsequent plots, time is scaled by the inverse of 
the mean-field interaction strength and hence is a dimension- 
less quantity; K3 is scaled by 1/vTV. The dashed lines give 
the exact results (Eqs. |12l) and the solid line gives the ap- 
proximate result (Eq. ( 13 i), which is accurate for large TV or 
small mean-field interaction time N\t- The inset shows K3 
for TV = 1000 in more detail. 



The fourth-order cumulant k 4 is plotted in Fig. [2] 
Again, for large photon numbers the cumulant ap- 
proaches a limiting scaling behaviour: 



which gives the same relative growth of 



K 4 



TV, 



(15) 



(16) 



although if the time is adjusted as a function of TV to keep 
the Kerr squeezing constant, the fourth-order cumulant 
decreases in proportion to the particle number. 

Figures [3] and [4] show the cumulants as a function of 
TV, for a fixed value of x^t — 25. One can clearly see 
two different regimes of behaviour, with the crossover 
between the two occurring just above TV ~ 10 4 . For « 4 
the number-state result is also plotted for comparison. 
One can see that K4/(Nt) 4 scales as described above for 
large N, but for small TV is limited to values of the order 
of corresponding number-state results (increasing with TV 
quadratically). This result suggests that a Kerr-squeezed 
state can be as non-Gaussian by this measure as the num- 
ber state, for sufficiently long interaction time. 
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FIG. 2: (colour online) The fourth-order cumulant, K4 of Y = 
X(n/2) in a rotating frame as a function of time, for various 
photon numbers ranging from 100 to 10 6 , as labeled. Time is 
scaled by the mean-field interaction strength, and K3 is scaled 
by 1/TV. For large TV or small mean-field interaction time 
TV%t, the results approach the same limiting curve oc t 4 . The 
inset shows K4 for JV = 1000 in more detail, both for the Kerr- 
squeezed state (dashed line) and for the number state (dotted 
line). 
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FIG. 3: (colour online) Third-order cumulant /ts(?) of the 
Kerr-squeezed state at fixed %Nt = 25 as a function of TV. The 
behavi our at TV > 10 5 reveals the scaling behaviour described 
by Eq. ( 13 1; the behaviour at TV < 10 4 is due to the saturation 
effect seen in the inset to Fig. [l] 



QUADRATURE VARIANCES, ENTANGLEMENT 
AND EINSTEIN-PODOLSKY- ROSEN 
CORRELATIONS 

Although entangled states have already been predicted 
to be produced by the intracavity nonlinear coupler [161 
117) . the linearisation process used to obtain the spectra 
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FIG. 4: (colour online) Fourth-order cumulant ^(S-) of the 
Kerr-squeezed state at fixed xXt = 25 as a function of TV 
(crosses). The squares give the fourth-order cumulant of a 
number state with the same number of photons. For TV < 10 3 
the Kerr-squeezed result saturates to that of the number state 

> 10 4 , the results follow the 



(see inset to Fig. [2]). For TV 
scaling as given in Eq. 
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in those cases forces Gaussian statistics on the outputs. 
Here, we wish to produce entangled states that maintain 
non-Gaussian statistics, so we will proceed by mixing the 
outputs of two Kerr oscillators on a beamsplitter [18] , as 
experimentally demonstrated by [TU]. We will now show 
the quadrature variances, as we need squeezed states in 
order to obtain entangled modes in the outputs. We note 
that, while this could be done by mixing one squeezed 
mode with vacuum, better results in terms of the degree 
of violation of the relevant inequalities are obtained by 
mixing two squeezed states. The Heisenberg uncertainty 
principle demands that 



V X 



> 1, 



(17) 



so that any quadrature with variance below one is 
squeezed. Fig. [5] shows the variances of the canoni- 
cal X = X(0) and Y = A(|) quadratures with time 
again scaled by N\t so that, apart from small- number 
effects, the same level of squeezing is obtained for differ- 
ent photon numbers. As for the cumulants, the results for 
N > 10 4 cannot be distinguished on this time scale. In 
fact, above N = 1000, the small-number effects only ap- 
pear as small differences amount of squeezing, so for the 
remainder of the paper, we quote results for N = 1000. 

Considering a beamsplitter with reflectivity 77 and la- 
belling the inputs by hi and a 2 and the outputs by b\ 
and b 2l we find 



h = - »?ai + \fr\a 2 - 



FIG. 5: (colour online) The variances in the X (solid line) 
and Y (dashed line) quadratures as a function of mean-field 
interaction time N\t- On this scale, the squeezing results for 
different numbers are indistinguishable for TV > 10 4 . 



For notational convenience, we will now make the sim- 
plification Xj, . Xj , and similarly for Yf, . . This allows 
us to define the variances of the beamsplitter outputs as, 



V(X 1 ) = r,V(X ai ) + (l-r,)V(Y a2 ), 

V(X 2 ) = {l-rj)V(Y ai )+nV{X a2 ), 

V&) = vViYaJ + il-^ViXa,), 

V(Y 2 ) - {l~r,)V{X ai )+r,V{Y a2 ). 



(19) 



Along with the covariances, 



v(x u x 2 ) = -y/^i-n) 



V(X ai ,Y ai ) + V(X a2 ,Y a2 ) 

V{Y U Y 2 ) = v^(l - V) [v(X ai ,Y ai ) + V(X a2 ,Y a2 ) 

(20) 

we now have all the expressions needed to calculate 
the quantities necessary to check for violation of the 
continuous- variable Duan-Simon inequality 19, 2fJj. For 
the purposes of this article, we define this as 



V(X 1 ±X 2 ) + V(Y 1 tY 2 ) >4, 



(21) 



(18) 



with any violation of this inequality being sufficient to 
demonstrate the presence of entanglement for a non- 
Gaussian state. The result for this correlation, with 
7] = 0.5, is shown in Fig. [6j The (red) dotted line gives 
the maximum violation, optimised over quadrature angle 
6. Clearly the outputs from two Kerr oscillators mixed 
on a beamsplitter can give a continuous-variable non- 
Gaussian entangled bipartite resource. 

As shown by Wiseman et al. [H] and Cavalcanti et 
al. |22) . the inseparability of the system density matrix 
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FIG. 6: (colour online) Duan-Simon correlations (LHS of 
Eq. ( |21[ )) after mixing on a 50 : 50 beamsplitter, as a function 
of interaction time x^i f° r N=1000. A value below 4 signifies 
entanglement. The (blue) solid line uses the canonical X\ — X2 
and Y1+Y2 quadratures and the green (dash-dotted) line uses 
the canonical X1+X2 and Y1 — Y2 quadratures, while the (red) 
dotted line is optimised for quadrature angle at each time. 



describes a set of states which includes within it subsets 
which are more deeply non-classical than evidenced by 
entanglement alone, such as those which demonstrate the 
Einstein-Podolsky- Rosen (EPR) paradox [23]. For our 
purposes here, we will use the inequality developed by 
Reid [24 , written as 



V inf {X bj )V inf (Y bj ) > 1, 



(22) 



where j = 1,2 and 



Vinf(X bj ) = V(X b] ) 



V inf (Y b .) = V(Y b .)- 



V(X b .,X bk ) 



V(X bk ) 

V(Y bp Y bk ) 
V(Y bk ) 



(23) 



From the expressions given above for the Duan-Simon 
criterion, it can be seen that all the moments necessary 
to calculate these expressions are available analytically. 
As shown in Fig. [7J the two modes after the beamsplitter 
exhibit a strong violation of the Reid inequality. 

Finally, we consider the skewness of the final entangled 
state. For a 50 : 50 beam splitter, the third- and fourth- 
order cumulants of the A"-quadrature in output port 1 



FIG. 7: (colour online) Reid EPR correlation (LHS of 
Eq. \22\ ) after mixing on a 50 : 50 beamsplitter as a func- 
tion of interaction time x^i f° r N = 1000. A value below 1 
signifies a demonstration of the EPR paradox. The dashed, 
lower line is optimised for quadrature angle, while the upper, 
solid line is for the canonical quadratures. 



can be shown to be: 

« 3 (*i) - -4 M*«i) - K 3 (Y a2 )) 

VO 

K4(*l) = ^(K4(^ a i) + K 4(^2)) + 4(X 1 ) K3 (^i), 



(24) 



which confirms that skewed inputs to a beamsplitter lead 
to skewed outputs, with cumulants generally of the same 
order of magnitude. 



CONCLUSIONS 

In summary, we have employed a simple model of 
nonlinear process to demonstrate that violations of the 
Duan-Simon and Reid entanglement criteria occur at the 
same time as significant departures from Gaussian be- 
haviour. For sufficiently long interaction time N\t, the 
nonlinear interaction will skew the distribution of the 
quadrature variables, leading to large third and fourth 
order cumulants. 

Moreover, such non-Gaussian entanglement occurs in 
regimes accessible to optical fibre experiments [5S] . How- 
ever, for accurate quantitative predictions, one would 
need to go beyond the single-mode model, to include the 
effects of pulse dynamics and extra noise sources, using 
the simulations methods, for example, that were used in 
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As with any coherent scheme employing the x^ 3 "* non- 
linearity of optical fibres, the intrinsic weakness of the 
nonlinearity itself can be a limiting factor. Besides using 
large photon-numbers, this factor may be overcome by 
use of electromagnetically induced transparency to pro- 
duce giant cross-Kerr nonlinear phase shifts and hence 
highly non-Gaussian quantum states |27j . On the other- 
hand, the presence of large numbers of photons makes 
optical fibre a very bright source of non-Gaussian entan- 
glement, which may well be a practical advantage over 
number-state schemes. 
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